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We study the effects of external offset charges on the phase diagram of Josephson junction arrays. 
Using the path integral approach, we provide a pedagogical derivation of the equation for the phase 
boundary line between the insulating and the superconducting phase within the mean-field theory 
approximation. For a uniform offset charge q — e the superconducting phase increases with respect 
to q = and a characteristic lobe structure appears in the phase diagram when the critical line is 
CI plotted as a function of q at fixed temperature. We review our analysis of the physically relevant 

situation where a Josephson network feels the effect of random offset charges. We observe that 
(— i | the Mott-insulating lobe structure of the phase diagram disappears for large variance (a > e) of 

Q ■ the offset charges probability distribution; with nearest-neighbor interactions, the insulating lobe 

O ' around q = e is destroyed even for small values of a. Finally, we study the case of random self- 

capacitances: here we observe that, until the variance of the distribution reaches a critical value, 
fS |i the superconducting phase increases in comparison to the situation in which all self-capacitances 

^ ' are equal. 
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Josepbson junction arrays (JJA) arc an ideal model to study a variety of phenomena such as phase transitions, 
frustration effects, and vortex dynamics [1]. The high level of accuracy reached in their experimental realization makes 
JJA relevant for different applications, ranging from possible implementations of quantum computation schemes [2] to 
, the realization of topological order [3]. Furthermore, particular geometries of JJA may induce novel classical coherent 
OO ' states [4]. 

^sO , JJA and granular superconductors, namely systems of metallic grains embedded in an insulator, become super- 
conducting in two steps [5]. First, at the bulk critical temperature, each grain develops a superconducting gap but 
the phases of the order parameter on different grains are uncorrelated. Then, at a lower temperature T c , the Cooper 
pair tunneling between grains gives rise to a long-range phase coherence and the system as a whole exhibits a phase 
transition to a superconducting state. The latter transition is governed by the competition between the Josephson 
tunneling, characterized by a Josephson coupling energy Ej, and the Coulomb interaction between Cooper pairs, 
described by a charging energy Ec- In classical junction arrays the Josephson coupling Ej is dominant and the 
transition separates a superconducting low temperature phase from a normal high temperature phase. When Ec is 
comparable to Ej (small grains) , charging effects give rise to a quantum dynamics and the energy cost of Cooper pair 
tunneling may be higher then the energy gained by the formation of a phase-coherent state. 

It is relevant to analyze in detail the effect of a background of external charges on the superconductor-insulator 
transition of a quantum JJA. Offset charges arise in real physical systems as a result of charged impurities or by 
application of a gate voltage between the array and the ground. In the former case, offset charges are distributed 
randomly on the lattice while in the latter case they play the role of a sort of chemical potential and their distribution 
can be uniform. Thus, offset charges may be regarded as effective charges qi, located at the sites of the lattice. When 
j_j • qi 7^ 2e, the offset charges cannot be eliminated by Cooper pair tunneling. As it will be clear, offset charges frustrate 
the attempts of the system to minimize the energy of the charge distribution of the ground state (for this reason they 
are also called frustration charges) . A large number of studies has by now been devoted to the analysis of the effects 
induced by offset charges on the (zero-temperature) quantum phase transition [6-9] and on the phase transition at 
finite T [10-12]. 
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The Hamiltonian commonly used to describe the Cooper pair tunneling in superconducting quantum networks 
defines the so called quantum phase model (QPM). In its most general form it is given by 

H = \ E(^ + '(Qj + «j) - e j E cos (^ - 

where is the phase of the superconducting order parameter at the grain i. Its conjugate variable m ([fi, ni] = z <5ij ) 
describes the number of Cooper pairs in the i-th superconducting grain. The symbol (ij) indicates a sum over nearest- 
neighbor grains only. The first term in the Hamiltonian (1) determines the electrostatic coupling between the Cooper 
pairs: Qi is the excess of charge due to Cooper pairs {Q\ = 2eri\) on the site i and Cy is the capacitance matrix. 
The diagonal elements of the inverse matrix Cy provide the charging energy: Eq — e 2 C ii 1 /2 = e 2 /2Co, where Co 
is the self-capacitance. The second term describes the hopping of Cooper pairs between neighboring sites (Ej is the 
Josephson energy) . The parameter a = zEj / AEc (where z is the coordination number) governs the superconductor- 
insulator transition: for a <C 1 the Josephson network is in the Mott-insulator state, while for a > 1 is in the 
superconducting one. An external gate voltage Vi contributes to the energy via the offset charge q\ = J2j Cy • This 
external voltage can be either applied to the ground plane or, more realistically, it may be induced by charges trapped 
in the substrate. In the latter situation q\ is naturally a random variable: we shall review our study of this physically 
interesting situation in Section IV. 

As it is well known [1], the QPM (1) is equivalent to the boson Hubbard model (BHM) in the limit of large particle 
numbers per junction. The BHM describes soft-core bosons hopping on a lattice [13] and is defined as 

ij i {ij) 

Here, b\ (&i) is the creation (annihilation) operator for bosons and n\ = b\b\ is the number operator. By writing the &i's 
in terms of their amplitudes and phases, and by neglecting the deviations of the amplitudes from the average number 
of Cooper pairs (n), we are lead to the QPM (1). An exact mapping between the two models has been derived in Ref. 
[14]. The hopping term is associated with the Josephson tunneling (2(n)t — ► Ej) whereas C/y — > 4e 2 C i T 1 describes 
the Coulomb interactions between bosons. The chemical potential in the BHM plays a role analogous to the external 
charge in the QPM [fi — ► q{). Thus, a QPM with random offset charges corresponds to a BHM with random on-site 
energies. 

In this paper we shall use a pedagogical approach to review the results obtained in our analysis [12,15] of the 
finite temperature superconductor-insulator transition in JJA. Using the finite temperature path integral approach, 
we provide an explicit derivation of the equation for the phase boundary for quantum JJA with offset charges and 
general capacitance matrix within the mean-field theory approximation [12,16]. Offset charges dramatically change 
the phase boundary: as it may be easily seen from Figs. 1 and 2, already with diagonal capacitance matrices, and 
for a uniform offset charge q = e, the superconducting phase increases with respect to the unfrustrated case. A 
lobe structure appears in the phase diagram when, at fixed temperature, the critical line is plotted against the offset 
charge q. In Fig. 3 we plot the value of a at which the transition occurs as a function of q for two different critical 
temperatures T c . As it is shown in Fig. 3, the critical line is 2e-periodic, and has a minimum in q = e, corresponding 
to the fact that for this value of q the superconducting region is maximal. 

Within the same mean-field approach, we shall also treat JJA with capacitive disorder [15], i.e., random offset 
charges and/or random self-capacitances. This is motivated by the fact that, in practical realizations of Josephson 
devices, one has to deal with capacitance disorder caused either by offset charge defects in the junctions or in the 
substrate (random offset charges) [17] or by imperfections in the construction of the devices, which may lead to 
random capacitances of the Josephson junctions. Although from a theoretical point of view charge and magnetic 
frustration are expected to be dual to each other, experimentally it is possible to tunc only the magnetic frustration 
in a controlled way. For this reason, it is widely believed that a challenging task for the theory is to develop reliable 
techniques to investigate the effects of random offset charges and capacitances on the phase structure of JJA. 

Random offset charges could be forced to vanish by using a gate for each superconducting island; however, this 
procedure works only for small networks, since in large arrays too many gate would be necessary, making impossible 
the cooling of the system at the desired temperatures. In Ref. [18] the case of uniform charge frustration was analyzed 
experimentally. It was observed a sensible variation (~ 40%) of the resistance between the unfrustrated and the fully 
frustrated array: however, it was impossible to quench the Coulomb blockade as it can be done in small arrays. 

During the last decade, much attention have been devoted to the Bose-Hubbard model with random on-site energies 
at T — [13,19-24]. As it has been already pointed out, the properties of the disordered BHM are closely related 
to those of JJA with random offset charges. Without disorder, the BHM exhibits two types of phases: a superfluid 
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phase and a Mott insulating phase. Strong disorder may lead to the existence of a third, intermediate gaplcss phase 
exhibiting an infinite supcrfhrid susceptibility and finite compressibility: the Bose glass (BG). It is by now known 
that, unless one introduces an ad hoc probability distribution, there is no BG in mean-field (MF) at T — [13]; the 
reason why - at T = - there is no BG phase in MF is better understood if one considers the infinite-range hopping 
limit of the Hamiltonian (2): tij — t/N for all the sites i ^ j, N being the total number of sites. At T = 0, when t = 
and (rii) = the system is, of course, in the insulating phase. If a small positive t is turned on, the infinite-ranged 
hopping allows the system to gain kinetic energy (with zero cost in the on-site potential energy) by moving bosons in 
unoccupied sites with the highest on-site energies; thus, the bosons are delocalized and the system becomes superfluid 
even for very small values of t. Of course, this argument does not hold at finite T; so far, there is no evidence for a 
Bose glass phase at finite temperature. 

In the following we shall review the results of Ref. [15], where the finite temperature phase diagrams of JJA with 
random offset charges and/or random self-capacitances have been obtained. 

The plan of the paper is the following: in Section II we use the coarse grained approach to compute the Ginzburg- 
Landau free energy for quantum JJA with charge frustration and a general Coulomb interaction matrix. The path 
integral providing the phase correlator needed to investigate the finite temperature critical behavior of the system is 
explicitly computed. 

In Section III, within the mean field theory approximation, we derive the analytical form of the critical line equation 
from the Ginzburg-Landau free energy. The phase diagram is drawn in the diagonal case for a generic uniform offset 
charge distribution. We then analyze the low temperature limit of a system with nearest neighbor interaction matrix 
when a uniform background of external charges q\ = e is considered. 

In Section IV we consider JJA with capacitive disorder (random offset charges and/or random self-capacitances): 
using the results of the previous Sections, we compute the free energy averaged over the disorder. Section IV. A is 
devoted to the study of the effects of random offset charges on the finite temperature phase diagram, with diagonal 
and nearest-neighbor capacitance matrices, while in Section IV. B we consider JJA with random self-capacitance. 

Finally, Section V is devoted to concluding remarks. 



II. PATH INTEGRAL APPROACH: THE GINZBURG-LANDAU FREE ENERGY 

The partition function for the frustrated model described by the Hamiltonian (1) is given by 

Z = Tre-^ = ^(^|e-^|^) (3) 

n 

where (3 = 1/ksT and the sum is extended only to states of charge 2e and thus with definite periodicity. 
In the functional approach Z reads 

z= |n^ cx p{-/ /3 ^(^w^w) } w 

where the Euclidean Lagrangian may be derived from 

1 = 5 © ? C "t f " (s) ? t « + E ' £ ~<" - *> (5) 

y i (ij) 

by replacing it/h — ► r. The path integral that one should compute is then given by: 

Z = jU^ ex P { r^[-i^C y || + i^ 9i | + f E^]} ^ 

i ij i ij 

where — oo < ipi < +oo, </?i(0) = tpi(0) + 2irni and 7y = 1 if i, j are nearest neighbors and equals zero otherwise (i.e., 
hopping term just between nearest neighbors). If 7y = 1 for all pairs i, j on the lattice, one is lead to the infinite-range 
hopping limit which provides a remarkable example of an exactly solvable MF theory [13]. The integers rii appearing 
in the boundary conditions take into account the 27r-periodicity of the states tjj n contributing to Eq. (3). 

In order to derive the Ginzburg-Landau free energy for the order parameter, it is convenient to carry out the 
integration over the phase variables by means of the Hubbard-Stratonovich procedure [25]: using the identity 
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eJ+rj = detTp J -Q ^2^. e -^r-^-j+^j (?) 

i 

the partition function may be rewritten as 

Z = J J|J^i^re^ dT ^*S iJ ^7 i 7V j ) e -s E//W ] ; (8) 

where the effective action for the auxiliary Hubbard-Stratonovich field ipi, SEff[ip], is given by 



2^ J 2e2e ^ y 2e 



The Hubbard-Stratonovich field ipi may be regarded as the order parameter for the insulator-superconductor phase 
transition since it turns out to be proportional to (e 4Vi ), as it can be easily seen from the classical equations of motion. 
From Eq. (9), the Ginzburg-Landau free-energy may be derived by integrating out the phase field (pi. 

Since the superconductor-insulator phase transition is expected to be second order [26], it may be safely assumed 
that, close to the onset of superconductivity, the order parameter ipi is small. One may then expand the effective 
action up to the second order in ipi, getting 

S Eff [iP}=S Eff [0}+ f dr [ P dT'G rs (T,T')A(r)r s (r') + --- , (10) 
Jo Jo 



where G rs is the phase correlator 

-i r2n r / 1 

= ( e ^-M-^ s (r')^ (11) 



G rs (T,T ) = — 



if, t ip*=0 



P 2 6ip r (r)5ip s (r') 
The partition function (8) can then be written as [12] 

Z = J Y[DiPiDiP?e-P F M, (12) 

i 

where F[ip] is the Ginzburg-Landau free energy; due to Eq. (10), up to the second order in one has 

m = / dT I dT ' ^ ^ (r)[7ij 1<5(r " T>) Gij(r ' t ' )] ^ {tI) - (13) 

In order to compute the phase correlator G rs one should evaluate the expectation value in Eq. (12) by means of the 
path integral over the phase variables (pi (r) . In performing the integration needed for the explicit evaluation of Eq. 
(13), one should take into account that the field configurations satisfy 

tpi(fi) - <pi{0) = 2-KUi. (14) 

It turns out very convenient to untwist the boundary conditions by decomposing the phase field in terms of a periodic 
field (pi(r) and a term linear in r which takes into account the boundary conditions (14); namely, one sets 

27T 

<Pi{r) = <pi{r) + —riiT , (15) 

with <pi(p3) — (pi(0). Summing over all the phases <Pi{r) amounts then to integrate over the periodic field <pi and to sum 
over the integers ni. As a result [12], the phase correlator factorizes as the product of a topological term depending 
on the integers rii and a nontopological one; namely, 



G rs (r;r') = ^ >- 

J^ iex p|/ ^r(-IC y ^|)| 
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For the sake of clarity, we briefly outline the computation of the factors appearing in the right-hand side of the 
previous equation. Firstly, one should compute the path integral 



(17) 



/ D0 ie J "MT)-i0 s (T') exp< | 


Jo dT (~h c uii^)^ 


/ D<pi exp<j 


ylo <M 2^U2e 2e)j 


\ 



Fourier transforming 4>i{r) according to 



with < r < and uo m = ^S-m, the numerator of Eq. (17) becomes 



= T> E ^ ( 18 ) 



/ ( ^ 
H II ^n^tn ex P 1 - E E Ci i U l<t>i,n<t>in 

i n=l *• P ij n=l 

. oo . N 

^ E (^,«e— - 4>l, n z- WnT ') + ^ (<Ar,o - 0.,o) + ccj . (19) 



n=l 

Upon integrating over the components </> r ,o, s o one gets a factor (5 



/oo \ / />oo />oo \ 

#i,o ( / #r,o / # s ,oe^-°^' o) -^-K 
.... -oo / \J—oo J —oo / 



(20) 



where X is an irrelevant divergent constant which cancels against the denominator. Using Eq. (20), Eq. (19) becomes 

oo „oc f ^ 

KS rs Yl / 11 Whndtpln exp I - -— E Cij w n0i,n0jj,„ 
n=l J -°° i [ P ij 

+ E <^M^ r - e-- T ') - E 0iJ^(e-^' e — ^')| . 

The multiple Gaussian integral is easily computed to give [12], up to an irrelevant constant which cancels against the 
denominator, 

oo „oo , . 

S *> H IT / d ^d^n exp E e "" T '> 

Ae 2 f3C u \ i . , ) 

i-^f-)/Me— * )} = 



i' Se 2 ^,. 1 ^ 1 - cosw„(r - r') 
^sexp]— _£( - ) 

n=l n 



where — /3 <t~t'<(3. By using the identity 



x v-^ / 4 4coso;„x\ „ . . 

N- 7 = Efe-— ) -0<x<0, ,2!) 

the first term of the right-hand side of Eq. (16) is found to be 

S rs exp|-2e 2 C rr 1 (jr - r'| - }■ (22) 

The sum over the integers in the topological factor in Eq. (16) is carried out by means of the Poisson resummation 
formula [27] 

| det G|i e~*'(™~ a ' iGij ( ,l ~ a ' J = ^ e -vm i (G- 1 ) i jm i -2wim i a i _ ^3) 
[ni] [mj] 

Due to Eq. (23), Eq. (16) becomes 

G rs (r,r') = <5 rs e- 2e2 ^ 1 l— 'I • 
y ,-E ii 2e2 ' 3C 'u 1 (™'+ft)("j+^)-Ei 4e2c H 1 K+li)(^-^) 

Mm] e 



y P E i j 2 / 3e2c u 1 ("'+li)("J+ft) 

Mm] e 

with ni taking all integer values and E[n-] being a sum over all the configurations. 
By means of a Euclidean-time Fourier transform, the fields tpi are defined as 



(24) 



V>i(r) = ^ViKK w " r 
n 

where w M are the Matsubara frequencies. As a consequence, the phase correlator Gy may be written as 



From Eq. (24) one can easily show that G rs (oj^;uj'^) is diagonal in the Matsubara frequencies and is given by 

Grs^jo^) = G r (w M ) • 5 rs ■ (5(0^ + uy) (26) 

with 

-IV £SL(„ I+ a)( n , + |l) 
1 — P y Z-(ij ^oo 

G r (w„) = V . (27) 

In Eq. (27) is given by 

z = £ e -e£ U 7^(»'+&)(«j+£) 

[m] 

with C/y = 4e 2 G i T 1 , i? c = e 2 G rr 1 /2 and y = k B T c /E C - In terms of Matsubara frequencies the Ginzburg-Landau free 
energy (13) becomes [28,11,12] 



1>lM ■ (28) 



Equation (28) is the pertinent starting point for the analysis of the phase boundary between the insulating and the 
superconducting phases in JJA with arbitrary capacitance matrix and with a generic charge frustration. 



6 



III. MEAN FIELD RESULTS 



In the following we shall illustrate the steps involved in the derivation [12] of the equation determining the phase 
boundary in the plane (a, KbT c / Eq), in mean-field theory and for a system with arbitrary capacitance matrix and 
a uniform distribution of offset charges. For this purpose, it is convenient to expand the fields ip-^oj^) and Gi(w M ) in 
terms of the vectors of the reciprocal lattice k. One has 



Gi(u 



e 



(29) 
(30) 



Moreover 



7ii 



(31) 



where 7 k is the inverse of the Fourier transform of the Josephson coupling strength 7y . As a consequence 

- 1 



E„e-' tp ' 



where p is a vector connecting two nearest neighbors sites. Expanding in k, one gets 



7 k " 



+ 



+ 



(32) 



where a is the lattice spacing. The first term in Eq. (32) provides the mean field theory approximation which, as 
expected, is exact in the limit of large coordination number. 
The Ginzburg-Landau free energy (28) reads 



/j-kk' 



2 _! Gk-k'(^u) 
^7 k ^kk ^ 



(33) 



Using Eq. (32) and keeping only terms of zero-th order in and k, one obtains the mean field theory approximation 
to the coefficient of the quadratic term of F: 



\M^)\ 2 ■ 



The equation for the phase boundary line then reads as 



l = z^G o (0), 



with 



(34) 



(35) 



(36) 



As evidenced in Ref. [12], the MF theory approximation amounts to neglect all the higher order terms in Eq. (33). 
Equation (35) determines the relation between T c and a at the phase boundary within the MF approximation. 

For a uniform distribution of offset charges Eq.(35) simplifies further since in Eq. (36) G r docs not depend on r. 
As a consequence, the phase boundary equation becomes 



1=a 'h ^c^^+dp 



(37) 
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with 

zEj 

a = TE c 

and 

In the following we shall illustrate through examples the physical implications of Eq. (37). 



A. Self Charging Model 

For a diagonal capacitance matrix, C7y = (Jyt/oo, one singles out only the self-interaction of plaquettes. Eq. (37) 
becomes 

- = <?(<?, 2/) (38) 
a 



where 



E, 



e 



We observe that Eq. (38) is invariant under g — > g + 2e and symmetric around q = n + e, where n is an integer: 

g(q + 2ne,y) = g(q,y); g [{2ne + e) + q, y] = g [(2ne + e) - q, y] . (40) 

In Fig. 1 we plot T c as a function of a. One sees that there is no superconductivity for a < 1. Due to the periodicity 
of Eq. (38) this holds for any integer q. For q = e, one gets 

e y + En=l e 



a = (41) 

4 +? y c -l , V +oo 1 -A (w+ i)a- ^ ^ 

4y e + Z^n=l l-4(n+i)2 e 



From Fig. 2, in which we plot T c vs. a, one sees that superconductivity is attained for all the values of a, since the 
superconducting order parameter at zero temperature is different from zero: a uniform offset charge q — e always 
favors superconductivity. In Fig. 3 we plot the phase boundary line as a function of the charge frustration, for a finite 
value of the critical temperature T c and also for T c — > 0. We observe that the limit T c — > is singular; nevertheless, 
taking this limit in the finite temperature Eq. (35), one finds an expression for the phase boundary line in good 
agreement with the results reported in Rcfs. [10,11,29]. 



B. Nondiagonal Capacitance Matrices 

In this Section we analyze the situation arising when the diagonal interaction matrix element Uqq and the nearest 
neighbor interaction matrix element Uo p = 9 Uqq arc nonzero. To do this, one should expand the critical line Eq. 
(37) for q = and small critical temperatures [30]: 

a=l+[|+2*(l- r ^)]e-5+... . 

When 9 > l/\/4 + 3z, the coefficient of the exponential e~ 4 / y is negative and the phase boundary line a — a{T c ) 
first bends to the left; when the critical temperature is high enough, it bends to the right, favoring the insulating 
phase. This is an indication for the possibility of observing reentrant superconductivity in these systems: fixing the 
parameter of the JJA (i.e., fixing a), by decreasing the temperature it is possible to go down from an insulating state 
to a superconducting one and then, further decreasing T c , to go back to the insulator. 
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As evidenced by Fishman and Stroud [30], the regime of physical interest is 6 < l/z; namely, when the capacitance 
matrix is invertible. Therefore, in dimensions D > 2, reentrance shall not occur with q = and nearest neighbor 
interaction matrix [30]. With q = e and a nearest neighbor interaction matrix, the situation is different [12]. We 
define 

*w = E|>+>* + 5> (42) 

i j 

the energy of a generic charge distribution on the lattice {rii}. Denoting with n? and n\ the charge distributions of 
the two lowest lying energy states and with E° and E 1 the corresponding energies, the low temperature expansion of 
Eq. (37) yields 

a- +E '" 1 " % <«) 

Eki + E "" 1 '-*[£i<-] + i)]' + " 

Independently on the explicit form of £/y, £7[ni] (for a square lattice in D dimensions) reaches its minimum value 
when (n? + \) = ±|(— l)* 1+,2 +— + tr > with ij (j = 1, ...,-D) the components of the lattice position vector i in units of 
the lattice spacing. 

For models with nearest-neighbor interaction, i.e., Uy/f/oo = <^ij + ^Sp^i+pJ with J2 P denoting summation over 
nearest neighbors, the first excited state has an energy E[n\] = E[n?] + z6, where S[n?], the ground state energy, is 
given by - z6)/A. 

With the above values of E[n9] and E[n\] and keeping only the leading order term in T c , Eq. (43) becomes 



where 



a = (1 - (1 - z&f) • (1 + a ie -y zf> + • • •) (44) 

/ l-(l-z9) 2 \ / z6) 2 \ 

Reentrant behavior at low temperature occurs when the coefficient of the exponential is negative, namely when a\ < 0. 
In Fig. 4 we plot T c versus a for 8 = 0.05 and z = 6. The resulting diagram exhibits reentrance in the insulating 
phase even for models with nearest neighbors interaction. For a detailed study of the lobe diagrams in JJA with 
non-diagonal capacitance matrices, see [1]. 



IV. CAPACITIVE DISORDER 



In this Section, we shall determine the finite temperature phase diagram of JJA with capacitive disorder (i.e., with 
random offset charges and/or random self-capacitances). To derive the phase boundary between the insulating and 
the superconducting phase, we shall use the path-integral approach for quantum JJA with offset charges and general 
capacitance matrices reviewed in the previous Sections. We find that charge disorder supports superconductivity and 
that the relative variations of the insulating and superconducting regions depend on the mean value q of the charge 
probability distribution: when q — 0, increasing the disorder leads to an enlargement of the superconducting phase. If 
the charge disorder is sufficiently strong (a > e), the lobe structure [1] disappears: in other words, the phase boundary 
line (and the correlation functions) do not depend any longer on q. In the following, we shall provide a quantitative 
analysis of this phenomenon. Also the randomness of the self-capacitances leads to remarkable effects, namely, the 
superconducting phase increases with respect to the case where disorder is not present. 

We shall consider several probability distributions which we expect to provide a realistic description of experimental 
situations. The low temperature behavior obtained by a pertinent extrapolation of our finite T results is consistent 
with the phase diagram obtained in Ref. [13]. 

For a given realization of the disorder, the Ginzburg-Landau free energy (i.e., the free energy near the transition) is 
given by Eq. (34) . We shall perform a quenched average, in which each of the random variables takes a unique value 
as the statistical variables fluctuate. This corresponds to taking the average of the logarithm of the partition function, 
i.e., the free energy The average of the free energy over all the possible realization of the disorder allows for the 
evaluation of the effect of a random charge frustration {<7j} or a random diagonal charging energy terms Ua = 4e 2 C ; T 1 . 
The pertinent starting point for the analysis of these situations is then 
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m = J d{x}p({x})FM 



(46) 



where P({X}) is a given probability distribution and d{X}P({X}) = JJ ; dqiP{q{) if one analyzes the effect of random 
offset charges or d{X}P({X}) = dUuP(Uu) for random charging energies. The random variables on different 
sites are taken to be independent. The phase boundary line between the insulating and the superconducting phase is 
determined by requiring that F = 0, which in turn leads to [15] 

l = z^-G . (47) 



A. Random offset charges 

In the following, we shall consider three different random offset charges probability distribution with mean q and 
width a . That is, a Gaussian distribution P(qi) = const ■ e~( q ''~ q ^ l 2a , a uniform distribution P(qi) = const between 
q — a and q + a and otherwise, and a sum of 5-like distributions P(qi) = J2 n Pn${Qi — ne ), with J2 n Pn = 1- For a 
diagonal capacitance matrix, Eq. (47) leads to 

-= f dqP(q)g(q,y) (48) 
a J 

with a = zEj/AE c , y = kBT c /E c and g(q, y) given by Eq. (39). If one considers the infinite-range hopping limit, one 
still gets Eq. (48) [15]. 

The results obtained from Eq. (48) with a Gaussian distribution are displayed in Fig. 5. One observes that, when 
q = 0, increasing a favors the superconducting phase while, when q = e, increasing a leads to the increase of the 
insulating phase. For large a (i.e. a > e), the phase boundary line is the same for all the values of q (in Fig. 5 the 
large a behavior is represented by the bold line). This is expected since, when a is large, the average free energy F 
does not depend any longer on q. 

As seen in the previous Sections an useful representation of the phase diagram is provided if one plots, at fixed 
critical temperature, the phase boundary line on the plane q — a. Without disorder one observes the lobe structure 
as discussed in Section III. In the presence of weak disorder and for T c — > 0, the lobes shrink, evidencing a decrease 
of the insulating phase: for a Gaussian (or unbounded) distribution the insulating phase completely disappears even 
for an arbitrarily weak disorder [13]. 

In Fig. 6 we plot the phase boundary line on the plane q — a at finite T c for the Gaussian and uniform distributions 
[15]. When the disorder increases the lobes flatten even at finite temperature and the same lobe structure is obtained 
from both distributions. In the limit T c — > 0, one recovers the result of Rcf. [13]. This can be easily seen if one 
observes that, at very low temperatures, for \q\ < e, one has from Eq. (39) 

9(q,y^0)= 4(^)2 - (49) 

Without disorder (a = 0), Eq. (48) simply gives a = 1 — A{q/2e) 2 . With the Gaussian distribution, since g has a pole 
in the half-integer value of the Cooper charge, the integral in Eq. (48) diverges and a — > 0, i.e., the lobes disappear 
for every value of a. As evidenced in Fig. 7, for a uniform distribution, when a > e, then a — > 0; when a < e, a — > 
only for e — a < q < e + a in agreement with [13]. 
Another interesting situation arises if one considers 

P(Qi) ne), (50) 

n 

with J2 n P n = 1- T ms corresponds to a random distribution of charges which are integer multiples of e and, actually, 
this is the most realistic situation for a random distribution. In fact, the probability distributions employed before 
should be viewed as fictitious continuous distributions, i.e., the properties of the overall distribution of charges (mean 
value and width) can be well approximated with a continuous distribution P(q). Inserting the probability distribution 
(50) in Eq. (48) one has 

~ = / d Q^2PnS(q - ne)g(q,y) = ^p n g{ne,y) + p n g(ne,y), 

n odd even 
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where ^2 odd E even ) is a sum restricted to odd (even) integer. From Eq. (40), one has g(2ne,y) = g(0,y) and 
g((2n + l)e, y) = g(e, y), which leads to 

- =Pog{0,y) +p e g(e,y), (51) 
a 

where po = ^2 even p n (p e — J2oddP n ) 1S tne probability that the offset charge q is an even (odd) integer multiple of e. 
In Fig. 8 we plot the phase boundary line (51) for po = p e = 1/2. 

We discuss now the effects induced by nondiagonal capacitance matrices in the limit T c — > [15]. The phase diagram 
without disorder becomes richer [1]; for concreteness, we shall consider on-site and a weaker nearest-neighbor (NN) 
interaction, i.e., the inverse capacitance matrix is restricted to diagonal and NN terms. If one defines 9 as the ratio 
between NN and diagonal terms, one should restrict only to z9 < 1 in order to insure the invertibility of the capacitance 
matrix [30]. Without disorder, at very low temperatures an insulating lobe around q = e appears [11]: the width of 
this lobe is z9/(l + z9). Putting W = l + z6, Eq. (39) for \q/2e\ < 1/2W gives g(q,y 0) = 1/[1 - 4W 2 (q/2e) 2 }; for 
1/2W < q/2e < 1 - (1/2W) it becomes, 



(2W£-l)(2W£-3) (2W(£-1) + 1)(2W(£-1) + 3)J- v ; 

In presence of disorder, Eq. (48) for a uniform distribution gives a = for (1/2W) — a < q < (1/2W) + a and 
1 - (1/2W) - u < q < 1 - (1/2W) + a. Thus, the lobe width decreases as (z6 - 2aW)/W. One sees that for 
a = z9/2W the insulating lobe around q = e disappears. This phenomenon is evidenced in Fig. 9. 



g(q,y -> 0) = -- 



B. Random Self-Capacitances 



In many instances, it may happen that the network's parameters are not uniform across the whole array: despite 
recent advances in fabrication techniques, with the use of submicron lithography, variation of junction parameters 
associated to the shape of the islands can be also of 20% [1]. Thus, it is relevant in many practical situations to study 
JJA with randomly distributed self-capacitances: this corresponds to have a random diagonal charging energy [15,31]. 
In Ref. [31] the effects of disorder on the spectrum of elementary excitations at low temperatures have been studied. 

In this Section we shall study JJA at finite temperature with uniform charge frustration q and random self- 
capacitance Cii, ignoring non-diagonal contributions to the capacitance matrix. The diagonal charging energy terms 
Ua arc related to the self-capacitances Cu via Uu = 4e 2 C i 7 1 and are assumed to be independently distributed according 
to the probability distribution P(Ua) oc e -( u u- u o) 2 /2a 2 _ The average charging energy is defined as Eq = Uo/8. The 
diagonal electrostatic contribution to the energy Un needs to be positive. 

By averaging the free energy (46), the equation for the phase boundary becomes [15] 

- = r dU^P-g(U,y) (53) 
a Jo U 

where now a = zEj/AEq, y = UbTc/Eq, and U = Uu/Uq; the function g(U,y) is given by 

The results of Eq. (53) are summarized in Figs. 10 and 11: when a is small, the superconducting phase increases 
in comparison to the situation in which all self-capacitances are equal: this is due to the factor 1/U in Eq. (53), 
which makes larger the contribution of junctions with charging energies less than Uo- An interesting observation is 
that, when q = e (maximum frustration induced by the external offset charges), the randomness does not modify 
considerably the phase diagram. This should be compared with the nonfrustrated case (q = 0), where randomness 
sensibly affects the phase diagram. 

The increase of the superconducting phase is due to a decrease of the effective value of the charging energy. This 
behavior occurs until a reaches a critical value (depending on the charge frustration and on the temperature), of 
order Uq- at this value of a the insulating region starts to increase. This is due to the asymmetry of the distribution, 
which has its peak in [7 , but only for positive values. This phenomenon is present also if one considers different 
distributions and is clearly seen in Fig. 12, where the phase diagram in the q-a plane for different values of the 
variance a is plotted. 
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V. CONCLUDING REMARKS 



In this paper, we reviewed the use of the path integral approach to finite temperature mean-field theory to analyze 
the effects induced by offset charges on the finite temperature phase diagram of Joscphson junction arrays. We 
provided, for a general Coulomb interaction matrix, the explicit derivation of the equation for the phase boundary 
line between the insulating and superconducting phase. 

The resulting phase diagram is drawn in the diagonal case for a generic uniform offset charge distribution q: with 
q = e, the superconducting phase increases with respect to q = 0, and the model exhibits superconductivity for 
all the values of a = zEj/AEc- An offset charge q = e tends to decrease the charging energy and thus favors the 
superconducting behavior even for small Joscphson energies. 

For the model with nearest neighbor inverse capacitance matrix and uniform offset charge q = e, we determined, 
in the low temperature expansion, the most relevant contributions to the equation for the phase boundary. For 
this purpose we explicitly constructed the charge distributions on the lattice corresponding to the lowest energies: a 
reentrant behavior is found even with a short ranged interaction. 

We also obtained the phase diagram at finite temperature of JJA with capacitive disorder. For a random distribution 
of offset charges with mean q and variance a, one has that for a > e, the phase boundary line coincides for any value 
of q and the lobe structure on the plane q — a disappears (a is the ratio between the Josephson and charging energies) . 
For very low temperatures there is agreement with the result of Ref. [13]. 

If one considers also a nearest- neighbor interaction, the insulating lobe around q = e, which arises in absence of 
disorder, is destroyed even for small values of a. For arrays with random charging energies, when the variance of 
the probability distribution is smaller than a critical value, the superconducting phase increases with respect to the 
situation in which all self-capacitances are equal. 

It is comforting to observe that the finite temperature MF theory approach developed in this paper provides results 
which arc in good agreement with those obtained by recent quantum Monte Carlo simulations [32] and by use of 
improved variational methods [33]. 
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FIG. 1. Phase diagram for the diagonal model without charge frustration. The critical temperature T c is in units of Jcb/Ec 
and q stands for the ratio zEj/AEc- The I and S indicate, respectively, insulating and superconducting phase. 
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FIG. 2. Phase diagram of the diagonal model with half-integer charge frustration q — e; T c is in units of ks/Ec- 
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FIG. 3. Lobe diagram for T c — > (solid line) and T c = ks/Ec (dashed line). As the Hamiltonian (1) is periodic in q with 
period 2e, the lobes are repeated on the q axis. 



a 



FIG. 4. Low temperature expansion of the critical line with a short-ranged inverse capacitance matrix. In the plot, z = 6 
and 6 = 0.05, where 9 is the ratio between nearest-neighbor and diagonal terms. 
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a 



FIG. 5. Critical line for random offset charges with Gaussian distribution (T c is in units of ks/Ec). The bold line is for 
a = e: for a > e no significant deviations from this line are observed for all the values of q. To the left (right), we plot q = e 
(q = 0); in the plot a/2e = 0.05 (dashed lines) and 0.3 (dotted lines). 



(a) 



(b) 



a 



a 



q/2e 



q/2e 



FIG. 6. Phase diagram with diagonal capacitances and random offset charges with uniform (a) and Gaussian (b) distribution. 
Top (bottom) of the figures: k B T/Ec = 1.5(0.5). We plot the cases a/2e = 0.05 (dotted lines), 0.3 (dashed lines)and 0.5 (solid 
lines). For large a/2e the phase boundary line is flat and it is the same for both distributions. 
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FIG. 7. Phase diagram at very low critical temperatures for a diagonal inverse capacitance matrix and random offset charges 
with uniform distribution. We plot a/2e = 0.1 (dotted line), 0.25 (dashed line) and 0.40 (solid line). 




FIG. 8. Phase diagram for random offset charges with the probability distribution given by Eq. (50) and diagonal capacitance 
matrix. In the plot we take po = p e = 1/2. 



1G 



0.7 



q/2e 

FIG. 9. Disappearance of the insulating lobe around q = e for a short-ranged inverse capacitance matrix. The phase diagram 
is plotted for T c — > and with random offset charges uniformly distributed. In the plot z9 is equal to 0.1 while cr/2e is 
respectively (dotted line), 0.015 (dashed line) and 0.03 (solid line). For this value of zO, the lobe disappears at a/2e = 0.045. 




FIG. 10. Phase diagram in the T c — a plane for random diagonal capacitance with Gaussian distribution and without charge 
frustration (T c is in units of fcs/f/o) while a/2e is 0.1 (solid line), 1 (dashed line) and 5 (dotted line). For each line, the 
superconducting (insulating) phase lies on the right (left). 
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FIG. 11. Phase diagram in the T c — a plane for random diagonal capacitance with Gaussian distribution and uniform offset 
charge q = e (T c is in units of Ub/Uo) while a/2e is 0.1 (solid line), 1 (dashed line) and 5 (dotted line). 




q/2e 

FIG. 12. Phase diagram in the q — a plane for random diagonal capacitance with Gaussian distribution and uniform offset 
charge at ksT/Uo = 0.1. We plot a/2e — 0.1 (solid line), 1 (dotted line) and 5 (dashed line). 
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